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We consider a quantum description for a spherically symmetric gravitational collapse of a massless 
scalar field. The effective scenario from loop quantum gravity is applied to a homogeneous interior 
spacetime. Classical singularity that arises at the final stage of our collapsing system, is resolved and 
replaced by a quantum bounce. Our main purpose is to investigate the evolution of trapped surfaces 
during the collapse in semiclassical regime. We show that, in this regime, there exists a threshold 
scale bellow which no horizon can form as collapse evolves towards the bounce. By employing the 
matching conditions at the boundary shell, quantum effects are carried out to the exterior region, 
leading to an improved Vaidya geometry. In addition, the effective mass loss emerging in this model 
predicts an outward energy flux from the interior quantum geometry regime. 
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I. INTRODUCTION 

There are two main aspects related to the final state 
of gravitational collapse of a star. The first is the sin- 
gularity formation, in the sense that as the radius of the 
star vanishes the matter energy density diverges at its 
centre. The second is the evolution of the horizon dur- 
ing the collapse. In the latter, if the trapped surfaces 
form as collapse proceeds, then the final singularity will 
be covered by the horizon and a black hole can form. 
Otherwise, if such trapped surfaces do not form as the 
collapse evolves, the radial null geodesies emerging from 
the singularity can reach the distant observer and hence, 
singularity will be naked [IH_E]. 

It is believed that the singularity problem will be 
overcome in a quantum theory of gravity. Loop 
quantum gravity (LQG) [MZ] j as a non-perturbative and 
background independent approach of quantum gravity, 
provides a fruitful ground to investigate the removal of 
singularities [5] . Results from the symmetry reduction of 
LQG (known as loop quantum cosmology (LQC)) lead to 
the conclusion that the spacetime singularities is resolved 
in quantum gravity [5]. On the other hand, the status 
of the classical singularities that arise at the late-time 
stages of the gravitational collapse, was studied in LQG 
[TOUTS] . Within this context, different matter such as a 
standard scalar field [TOUT!?] or a tachyon field [T31 [TSJ 
have been considered as the collapsing matter source. 
Therein, by employing the quantum gravity effects (such 
as the inverse triad correction imported from LQG), it 
was shown that the geometry of spacetime near the clas- 
sical singularity is regular. Furthermore, some novel fea- 
tures, such as evaporation of horizons in the presence of 
quantum gravity effects was studied in Refs. [TUJ []__] . In 
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addition, it was shown in Ref. that (inverse triad 
modifications) quantum gravity effects predict a critical 
threshold scale for the horizon formation which may lead 
to the formation of very small non-singular astrophysical 
black holes. 

Improved LQC dynamics with a massless scalar field 
has been developed in the last few years [TT] [TH] (see also 
Ref. [H]). Concerning the physical implications of the 
singularity resolution in this model, it was shown that 
the classical big-bang singularity is resolved and replaced 
with a quantum bounce in LQC [171 118] . In view of 
these elements, it is expected that the singularity arises 
at the end state of the gravitational collapse would be 
also resolved and replaced by a bounce. However, the 
question further arises is how loop quantum effects can 
indeed affect on the emergence of trapped surfaces in this 
model. 

A trapped surface in classical general relativity is 
denned as a compact 2-dimensional smooth space-like 
submanifold of spacetime such that the families of out- 
going as well as ingoing future-pointing null normal 
geodesies are contracting [2J. When quantum geometry 
becomes relevant, geodesies stop and the classical state- 
ments, based on the properties of geodesies on a differ- 
ential geometry, are not valid anymore and has to be re- 
placed by something more appropriate in quantum the- 
ory. However, in order to extract physical information 
out of the theory it is of interest to employ the effective 
theory of LQG. Quantum dynamics of LQC, including 
a holonomy corrections, can be approximated by an ef- 
fective continuous equations of motion, which results in 
an effective theory of LQC 18 . The effective theory 
shares the form of the classical theory, but contains cor- 
rection terms from the quantum theory. Furthermore, in 
semiclassical regime, this scenario agrees with classical 
general relativity. Within such a scenario, the classical 
statements (such as singularities and trapped surfaces), 
looks promising to be developed in quantum theory. 

To that purpose, we apply the recent results of the 
effective theory of LQC in the resolution of the singular- 
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ities related to the gravitational collapse of a star. More 
precisely, we consider in this paper, a spherically sym- 
metric framework for gravitational collapse, whose mat- 
ter content includes a scalar field, and study the physical 
conditions for the trapped surface formation (see also 
Refs. p~0l - ll2| ) . Then we consider the effective dynam- 
ics scenario of LQC to study the collapsing system. We 
investigate whether such quantum gravitational correc- 
tions can affect the dynamical evolution of the collapse 
and determine the formation or avoidance of the trapped 
surfaces. 

The context of this paper is organised as follows. In 
section [TTJ we provide the background scenario, by in- 
troducing the choice of spacetime geometry for the col- 
lapsing system as well as the matter source. In partic- 
ular, we consider a flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) interior spacetime to be matched to the 
generalised Vaidya geometry at the boundary of matter. 
The matter source is considered to be a homogeneous 
and massless scalar field. Then, we describe the model in 
terms of suitable variables applicable in LQG. In section 
III[ we study the quantum gravitational collapse of our 
model by employing the effective dynamics scenario of 
LQC. In section [TV] we investigate how quantum effects 
influence the evolution of the trapped surfaces as the col- 
lapse evolves. Furthermore, by employing the matching 
conditions at the boundary of the collapsing cloud, we 
show that the quantum effects from the semiclassical in- 
terior region are carried out to the exterior region which 
leads to an effective exterior geometry. Depending on the 
dynamics of the trapped surfaces on the interior quantum 
geometry, we have scenarios where the final quantum 
bounce is visible to a distant observer or is covered by 
a non-singular quantum black hole horizon. Finally, we 
present the conclusion of our results in section |Vj 



II. GRAVITATIONAL COLLAPSE WITH A 
SCALAR FIELD 



The physical model we consider is that of a collapsing 
spherical body (star), whose matter is described by a 
homogeneous massless scalar field 1 . The matter is con- 
fined to a spherically symmetric region, whose coordin- 
ates is considered as (t,r,9,(f>). In order to describe the 
whole spacetime structure, we assume that the spacetime 
is splitted into two parts: The interior region which con- 
stitutes the matter field; and the exterior region which 
must be matched to the interior one at the boundary with 
the coordinate radius r = rj,. 



1 This belongs to the general type I matter, which also includes 
dust and perfect fluids. The main property is that the energy 
momentum tensor admits one time-like and three space-like ei- 
genvectors pQ: In the case at hand, they are the energy density 
p and the pressure p. 



The geometry in the interior region can be modeled by 
a metric of the FLRW family [10 1 H3 ] : 



g- v dx^dx v = -dt 2 + a 2 (t) (dr 2 + rW) , (2.1) 

where dft 2 is the standard line element on the unit two 
sphere. We can identify every shell with its coordinate 
radius r. The physical radius of such a shell is given by 



R{t,r) := a(t)r, 



(2.2) 



called the area radius. In light of the canonical analysis, 
which differentiates r and t, it is a good idea to fix the 
coordinate r and regard R(t, r) as a function on the grav- 
itational phase space. 

The Einstein's equations for the interior region can be 
presented [4] as 



8nG p = 



R 2 R, r 



8irG p = - 



F 
R^R 



R 2 = ~ , (2.3) 



where Go is the Newton constant. A 'dot' and ',r' de- 
note the differentiation with respect to the proper time t, 
and to the coordinate r, respectively. The mass function 
F(t, r) is the total gravitational mass within the shell 
labelled by r. For a massless scalar field <fi, the energy 
density p and the pressure p coincide, and can be ex- 
pressed in terms of the matter dynamical variables as 
p = p = 7T 2 /2a 6 . 

In order to investigate the geometry of trapped surfaces 
inside the star, it is convenient to study the behaviour 
of the radial null geodesies emerging from the interior 
spacetime. Let us introduce the null coordinates 



di ± = — t= [dt T a(t)dr] . 
V2 



(2.4) 



Then, the interior metric (2.1) can be recast into the 
double null form [2H1 as 



9p V 



dx^dx" 



-2d£ + dC +R 2 dn 2 . 



(2.5) 



The radial null geodesies are then obtained by solving 
g~ v dx^dx v = with the condition dfl 2 = 0. From here 
we deduce that there exist two kinds of null geodesies, 
corresponding to £ + = const., and £~ = const. We 
can compute the expansion parameters 9± for these 
geodesies, given by [2U] 



R 



R 



dtT 



W) 



R. 



(2.6) 



measuring whether the bundle of null rays normal to the 
sphere is diverging (9± > 0) or converging (9± < 0); R is 
the Ricci scalar curvature for the metric (2.1). Introdu- 
cing a new parameter 0(i, r) :— 9+9-, we are able to say 
whether the spacetime is respectively, trapped, untrapped 
or marginally trapped, depending on 



e(t,r)>0, e(i,r)<0, 9(i,r) = 0. 



(2.7) 
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The third case in Eq. (2.7) characterises the outermost 



boundary of the trapped region, namely, the 'apparent 
horizon" which corresponds to the equation R 2 — 1. Us- 
ing Eq. (2.6) we can write 0(t, r) as 



e 



1 R 2 



1 



2 \ R 2 R 2 



(2.8) 



This is also to be thought of as a function of the phase 
space, for every fixed shell r. Specifically, we will be most 
interested in the boundary shell, r = r;,, which bounds 
the support of matter. In that case, we define 



6 (f) :=&(t,r b ) = 



1 



1 



(2.9) 



Since we are mainly interested in the trapped surfaces 
eventual formation due to the gravitational collapse of 
the interior spacetime, we assume that the star is not 
trapped from the initial configuration at ti\ in other 
words, 6(i,, r) < for all shells < r < r^. 

To model the exterior geometry, we choose a metric 
of the Vaidya family 2 . Written in advanced Eddington- 
Finkelstein coordinates (v, r v , 8, </)), it has the form [21] 



g+ v dx»dx v 



(1 - 2Af(v)G /r v ) dv 2 
- 2dvdr v +r 2 dn 2 . 



(2-10) 



where M(v) is a generic function of v, which is fixed by 
matching the Eq. (2.10) with Eq. ( |2.1[ ) at the boundary 
r = rb (for a discussion, see [TU1 113]). The matching 
condition is given by matching the area radius at the 
boundary £ [4]: 

r v (v) = R{r h ,t)=r h a(t), (2.11) 
together with the first and second fundamental forms 



R,r + rbd 
1 _ £ 

1 R 



F{t,r h ) = 2M(r v ,v)G , 

pi 

G M(r v ,v) jT . v = — + rlaa . 



(2.12) 
(2.13) 
(2.14) 



It should be noted that the singularity formation at a = 
is independent of these matching conditions. 

Matching the exterior generalised Vaidya geometry to 
the interior spacetime plays two important roles in a col- 
lapsing process: Firstly, it allows the matter to be radi- 
ated away as collapse evolves; and secondly, it enables 
the study of formation and evolution of horizon during 
the collapse. The second aspect is particularly import- 
ant; indeed, formation of a black hole as the end state 



2 This is a generalization of Schwarzschild metric which accounts 
for possible matter emissions, and realizes the astrophysical real- 
istic case of a star surrounded by a radiating zone 12 1 1 . 



of a collapsing star indicates that there exists a moment 
when an apparent horizon develops inside the cloud, and 
thence all the matter content collapses inside such hori- 
zon. On the other hand, if a black hole is not the end 
state, it means that trapped surfaces never develop at any 
stage of the collapse, and hence apparent horizons never 
form inside the star. In this paper we are mainly con- 
cerned with whether or not trapped surfaces can form in 
the interior region, once quantum gravity corrections are 
taken into account (see section III). For this reason, we 



will consider only the quantum correction to the space- 
time inside the collapsing star, however, this quantum 
effect can be carried by the exterior geometry through 
the matching conditions. 



III. QUANTUM GRAVITATIONAL COLLAPSE 

In this section, we discuss the quantum geometry cor- 
rections to the classical evolution presented in previous 
section. To do this, we need to replace the phase space 
variables of the collapsing spacetime with Ashtekar- 
Barbero variables, more suitable for LQG quantization 
scheme. For the symmetry reduction of the interior 
FLRW type (2.1), the phase space of gravity, r grav , is 



two dimensional and coordinatised by the scale factor a 
and its conjugate momentum n a = — (6/87rGo)ad, which 
satisfies the Poisson algebra {a,ir a } — 1. To be able to 
apply LQC techniques, we make a canonical transform- 
ation to a more suitable pair of variables: Let us define 
the oriented volume as 



a 3 /a, 



(3.1) 



where a = 2^^\/~Ki\ x (with A = 4v / 37r7£| 1 being the 
so called area gap, and 7 is the Immirzi parameter) [5]. 
We easily find its conjugate momentum b = —j\/~Ka/a, 
with which it satisfies the Poisson algebra {v,b} = 2. 
Notice that, for our collapsing model d < 0, so that 
b > 0. In terms of these variables, the (interior) gravita- 
tional Hamiltonian constraint can be written as C grav = 
- (37rG /2a) \v\b 2 . 

As for the (massless scalar field) matter, the Hamilto- 
nian analysis for the homogeneous case, results in a two 
dimensional phase space of matter, r mat t, coordinatized 
by <j> and 71^ = av(f> satisfying {</>, n^} = 1. The matter 
Hamiltonian constraint is simply G ma tt = tt^/ (2a|u|). 
The total Hamiltonian constraint is then, given by 



C = C„ 



G n 



37rGo 
2a 



v\b 2 



2a\v\ 



(3.2) 



Notice that, since the scalar field cj> does not enter in 
the expression of the constraint equation (3.2), its mo- 
mentum 7T0 is a constant of motion. 

At this classical level, it is possible to solve the 
Hamilton equation for the field <p analytically: The equa- 
tion is 4> — {4>i C} — (a\v\). We can also write <j> 
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as <j) — <f>'v = 4>'{ V ,C} — — (j>'(6nGo/a)\v\b (where a 
prime denotes the derivation with respect to v), which 
using the constraint C = 0, can be rewritten as <p — 
(f/ir^y/Yl'KGo/a. Hence, we get 



1 



dv 



whose solution is 



\v\V12ttGo ' 



lnf- 
VTh^ \v 



(3.3) 



(3.4) 



Here (<p a , v ) are integration constants describing the ini- 
tial conditions for the collapsing star at the t = t space 
slice. We see that for (j> — > oo, it is v = 0, i.e., the volume 
of the interior region vanishes; but since the matter is 
contained in such region, the energy density of the cloud 
diverges, producing a physical singularity. We will see 
that once quantum corrections are taken into account, 
the situation changes completely. 



A. Loop quantization 

Following Dirac program for quantization of con- 
strained systems, we introduce a kinematical Hilbert 
space Hkin mirroring at the quantum level the symplectic 
structure of the phase space T = r grav x r mat t, and then 
restrict to the kernel of the operator representing the 



constraint C»0. Let us choose %kin = % g rav ® %matt, 
where 7i gI av = L 2 



c^Bohr)and"H = L 2 (R,d<^) 0. The 
main difference with respect to the canonical approach 
is the choice of a polymeric representation for the grav- 
itational sector, rather than the usual Schroedinger one. 
Indeed, the Hilbert space of gravity is the space of square 
integrable functions on the Bohr compactification of the 
real line, R, with respect to the Haar measure d/^Boiu- 
defined on it. To our purposes, we only need to know 
that on Hgrav the volume operator v is well defined, and 
its eigenstates form a complete orthonormal basis 3 : 



= <5„ 



(3.5) 



The result of using R instead of R is that there is no 
well defined operator representing b. The reason of this 
can be somehow found in the fundamental discreteness 
of space, which in turn implies that derivatives (such as 
b in the standard Schroedinger representation of the al- 
gebra {v,b} = 2) do not exist, as their definition would 
require infinitesimal displacements. Following this inter- 
pretation, we expect no problem for finite displacements, 
i.e., for the exponentiated version of b. Indeed, on "H grav 



3 Notice that the orthonormality is with respect to the Kronecker 
delta, 5 g , „/ . 



it is well defined the operator N\ := exp(iA6/2), where 
A is a constant. So that, 



N x \v) 



A). 



(3.6) 



As for the matter part, it is clear from our choice of Hil- 
bert space 'H ma tt that we use the standard Schroedinger 
representation, with <p acting as multiplicative operator 
and 7?0 acting as derivative operator. 

Once the action of fundamental operators is defined, 
any other phase space function / can be represented 
formally by 



f(v, 6, (/>, 7T ) = f(v, N, (j), 7?^). 



(3.7) 



Of course, one is still faced with ordering ambiguities. 
As an example, we construct the operator corresponding 
to 9. First, we rewrite Eq. (2.8) in terms of the new 
fundamental variables: 



0(t,r) 



1 



I 



7 2 A 



(3.8) 



We immediately meet two problems: (i) the function b 
appears, which has no correspondence at the quantum 
level; (ii) the function v~ 2 / 3 is non polynomial in v. In 
order to deal with the first problem, one observes that 
classically the following relation holds: 



1 

27 



iXb 



-i\b\ 



= Xb + 0(X 3 



(3.9) 



So, for small values of A, we may write b = (e 



Xb 



-iXb 



)/2i\. Now, the functions appearing on the right 
hand side of this expression have well defined quantum 
analogue, thus, we can write 



27A 



(3.10) 



To deal with the second problem, one can use the so- 
called Thiemann's trick [32]: Since classically v~ 2 ^ 3 — 
-(3i/\){v 1/3 ,e lXb/2 }e- iXb/2 , at the quantum level one 
has 



-2/3 



71/3 



A 



(3.11) 



Therefore, when we need the quantum correspondence 
of ii~ 2 / 3 in any expression, we can simply use the right 
hand side of the equality (3.11 ) instead. So, at the end of 



the day, the operator corresponding to O takes the form 



1 



4A 2 7 2 A 
3 

Aa 2 / 3 r 2 



Nx -2 



/3 - N~xv 1/3 N X 



(3.12) 



The most important operator to be implemented on Hkin 
is the Hamiltonian constraint (3.2 ). Indeed, the final step 



in Dirac quantization is to restrict our attention to the 
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states belonging to KerC, since they are the physical 
states. Looking at its classical form, Eq. (3.2), we again 



equations 



see the same problems that attacked 0. Despite the fact 
that we have a closed form for this operator, we will not 
be using it in this paper. Indeed, our propose herein this 
paper is to study the dynamics of the trapped surfaces 
evolving during the collapse by employing the classical 
statements 4 . To this aim, it is sufficient to study our sys- 
tem within the semiclassical regime, by resorting to the 
"effective dynamics" procedure |18j . In other words, we 
will forget about the fundamental quantum theory, and 
consider only the first order modifications to the classical 
theory, thereby obtaining 'effective Hamiltonian' equa- 
tions of motion [23 • 



B. Effective theory 

A well defined treatment exists, based on the exist- 
ence of a "minimal loop" around which the holonomy 
is computed, and used as the best approximation for b. 
However, for our purposes, it is enough to adopt the "ef- 
fective dynamics" methods. Since we are interested in 
an effective classical theory, we will not turn C into an 
operator, but rather consider the classical modifications 
that would mimic the quantum behaviour at first order. 
The only such a modification we will take into account is 
the so-called holonomy correction, which consists in re- 
placing b with the phase space function derived by the 
classical observation that 



I 

2iX 



AXb 



-iXb\ 



sin(A&) 



(3.13) 



obviously inspired by Eq. (3.10). Since this is the only 



change, it is expected that the classical theory is re- 
covered for small b, i.e., as long as the Hubble parameter, 
a I a is small; in other words, the difference in the dy- 
namics of observables will be non-negligible only in the 
vicinity of the classical singularity. 

we get 



Plugging the substitution (3.13) in Eq. 
the effective Hamiltonian [18, 23 : 



(3.2) 



C e ff = 



37rG 0, I • 2 / \ t\ i n - 

\v\ sin (Xb) + 



2aX 2 



2a\v\ 



(3.14) 



Now, the dynamics of the fundamental variables 
v, 6, 4>, 7T0 is obtained by solving the system of Hamilton 



4 The definitions of singularity and trapped surfaces are based 
on the evolution of the geodesic curves on spacetime manifolds. 
Therefore, in quantum theory these definitions are not valid any- 
more and can not be employed at the quantum level. 



v = {v,C e g} — — ° '|«| sin(A6) cos(A&), 



b = {b,c cS } 

4> = {4>, Ceff} = 



3ttGo . 2/ u \ , <f> 

sm (Xb) H hr 

ctv A 

fl> = 0. 



A 2 a 
a\v\ ' 



(3.15) 



These equations are not independent, since C e g = (i.e., 
the motions are curves in T lying on the constraint sur- 
face defined by imposing the vanishing of the effective 
Hamiltonian). 

Recalling that a — (va) 1 ^ 3 , one sees that H = a/ a — 
v/3v, so using the effective Hamilton equation for v 
in Eq. (3.15), and the constraint C c g = 0, given by 
Eq. ( |3. 14 1, one finds the effective Friedmann equation 



H 2 



8ttGq 



P 



(3.16) 



where p = 7r 2 /(2a 6 ) and p cr = 3/(87rG 7 2 A 2 ) ~ 0.41p P i. 
Eq. (3.16) implies that the classical energy density p is 
limited to the interval po < p < p C r having an upper 
bound at p CI . Notice that, po <C p cr is the energy density 
of the star at the initial configuration, t — 0, where po = 
*$/(2og). 

On an effective geometry, in order to compare with 
the classical general relativity, the modified Friedmann 
equation (3.16) is obtained by considering an effective 



Newton's constant given as G c g = G$ (1 — p/p cr ), where 
G is the low energy Newton's constant. In other words, 
it is interesting to consider the dynamical trajectories 
of the collapsing system to be given by the Hubble rate 
H = a 2 /a 2 , whereas the energy density is effective in a 
quantum geometry regime and is provided by the modi- 
fied Friedmann equation (3.16). Hence, we can write 



8irG 



efit 



8ttG 



3 r 3 
So, the effective energy density reads 

_P_ 

Pa 



PcS- 



Peff := p 1 - 



(3.17) 



(3.18) 



We see that quantum geometry effects lead to a modific- 
ation proportional to p 2 , which becomes important when 
the energy density becomes comparable to p CI . Further- 
more, in the limit p — > p CI , the Hubble rate vanishes; 
the classical singularity is thus replaced by a bounce (cf. 
figure [T]). In this sense, quantum geometry leads to a 
singularity resolution at the late-time evolution of the 
collapse. Notice that, in the limit p <§; p cr , the standard 
Friedmann equation is recovered. 

The classical Friedmann equation corresponds to the 
last relation in the classical Einstein's field equation ( 2.3 1 , 
given by H 2 = F/R 3 . Consequently, the effective theory 
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Figure 
(dashed 
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1: Behaviours of the scale factor a(t) in the classical 
curve) and effective (solid curve) cases. 



of LQC leads to a modification for the mass function in 
Eq. (2.3). Therefore, the mass function F(t,r) can be 



modified on the effective geometry as 

8irG 3 

^cff — g Peff-K • 



(3.19) 



The phase space trajectories are considered classical, 
whereas the matter content is assumed to be quantized. 
In the classical limit, as p c s —> p, the effective mass 



function reduces to the classical F 

F = (8TTG /3V2) P y bv fp. Eq 
p Q < p < p, 



lTrG /3)pR 3 or 



(3.191 shows that, as 
< p < p CI on the interior quantum geometry, the ef- 
fective mass function remains finite during the collapse. 
It is convenient to rewrite the Eq. (3.19) as 



F 



off 



F 1 



where F c , 



P< 



R 3 



So that, we can also write 
F 2 



F 2 
F 2 . 



— 8ttGo 2„3 
~ 3^2 * 



(3.20) 



is a constant. 



Pc 



(3.21) 



Eq. (3.20) shows that, mass function F changes in the 
interval Fq < F < F cr along with the collapse dynamical 
evolution. Nevertheless, if on the initial condition the 
mass function starts with F < R, then as the scale factor- 
decreases, the mass function also decreases and becomes 
zero at the bounce; consequently, there will be no trapped 
surfaces forming. This result is consistent with the case 
presented in the previous section. 



IV. EFFECTIVE DYNAMICS OF HORIZONS 
AND EXTERIOR GEOMETRY 

To discuss the dynamics of the trapped region in the 
perspective of the effective dynamics scenario, particular 
importance is played by the function 0b given by Eq. 



(2.9). Therein, by replacing H with the effective Fried- 



mann equation (3.16) we get 



AttGq 



1 



2a 2 ? 



(4.1) 



We can study the behaviour of the effective <db as a 
function of scale factor a. Since the cloud is initially 
untrapped, at the initial space slice (i.e., for large a) 
@b(t — 0) is negative. Then, as the collapse proceeds 
(i.e., a decreasing towards the collapse center), the second 
term always increases, whereas the first term (the Hubble 
rate) increases, whose effect being dominant until the en- 
ergy density reaches to p — 0.4p cr (cf. figure [2]); at this 
point 0fc reaches its maximum, max : 



©max = 0.327tGo/O c 



(4.2) 



>' b 



Figure [2] shows the behaviour of ©f, against the scale 
factor a for the different choices of the initial conditions. 
Solid curve represents the trajectories provided by the 
effective dynamics gravitational collapse. Dashed curve 
shows the classical trajectories which coincides the effect- 
ive equations for large values of a. Equation of apparent 
horizon on the effective geometry can be obtained by set- 
ting @b = 0. Therefore, depending on the initial condi- 
tions, in particular on the choice of the rj,, three cases can 
be evaluated. The left plot (solid curve) indicates a con- 
dition in which an untrapped interior spacetime without 
any horizon forming, whereas two others show trapped 
regions; the middle and the right figures correspond to 
one and two horizons formation, respectively. In addi- 
tion, denoting by dashed curves, always only one horizon 
can form classically. 

The scale factor a max , corresponding to max reads 

i 

/ -7T 2 \ 6 

(4.3) 



0.8/Jc 



Notice that this value is independent of r^, so it is the 
same for any shell. The minimum value of the scale factor 
is fixed by the requirement that the Hubble rate vanishes, 
i.e., p = p CI where the collapse hits a bounce as its final 



state. Using the relation p - 
putes the scale factor 



2p c 



7T?/(2(i ), one easily com- 



(2.5) 



1/6, 



(4.4) 



The modified Friedmann equation (3.16) allows to 
compute the maximum speed of the collapse |a| max , cor- 
responding to the scale factor a max : 



= 0.32ttG, 
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(4.5) 



On the other hand, by Eq. (2.8) we can determine the 
speed of the collapse, |d| ah , at which horizons do form, 



Figure 2: Behaviours of @b{a) in the classical (dashed curve) and quantum (solid curve) cases, for different values of r;, (from left to 
right: < r*, ri, = r* and > r*). Where r* defines the threshold radius. 



by simply setting 8 = 0. Then, we get R 2 
the boundary shell, gives 



a AH 



1 , which for 



(4.6) 



When the speed of collapse, |d|, reaches the value 1/rb, 
then an apparent horizon forms. Thus, if the maximum 
speed |d| max is lower than the critical speed |<z|aHj no 
horizon can form. Let us introduce a threshold radius r*, 
defined by 



1 



(4.7) 



We see that, if r b < r*, then no horizon can form at any 
stage of the collapse. The case r b — r* corresponds to the 
formation of a dynamical horizon at the boundary of the 
two spacetime regions Finally, for the case r b > r* 
two horizons will form, one inside and the other outside 
of the collapsing matter (cf. figure [3]) [TT] . 

The total mass m inside the collapsing star is given by 
the integral of the energy density on a sphere of radius 
R b = ar b : 



rn = a p<f,d x 



27T7T 



^ r 3 
l 



(4.8) 



Thus, we can translate the condition r b > r*, on the 
radius, into a condition on the mass of star. Let us define 
a m* to be the mass of a star with the shell radius jv 



2?r tt 



(4.9) 




Figure 3: The speed of collapse, |d|, with respect to the scale 
factor a, in classical (dashed curve) and semiclassical (solid curve) 
regimes. 



So far we have analysed the interior spacetime of the 
collapse in the presence of the quantum gravity effects. 
This quantum effects is expected to be carried out to 
the exterior geometry by using the matching conditions 
on the boundary r b of two regions. In the following we 
will focus on the main physical consequences that can 
emerge from this scenario in order to predict the possible 
exterior geometry for the collapse. To this aim, we will 
consider two cases: First, we will investigate the fate of 
the collapsing star whose mass is less than the threshold 
mass to*. Secondly, we will study the exterior geometry 
for the case in which the mass of the star is bigger than 
the threshold mass whereby a black hole can form. 



Therefore, in order for trapped surfaces to form, the 
mass, to, of a collapsing star must be bigger than the 
threshold mass to*. In other words, the minimum mass 
of the final black hole must be bigger than the threshold 
mass, to > to*. This inequality cannot be satisfied in 
the case of micro black holes, since in that case to is 
extremely small. Thus, for such black holes (which are 
expected to form in the primordial universe, as well as 
in LHC), trapped surfaces never form, and so informa- 
tion about the interior structure can travel outside, to 
the distant observer. 



A. Outward flux of energy 

For a collapsing star whose mass is less than to*, we 
study the mass loss using the quantum corrections to the 
interior spacetime due to the effective theory we discussed 
before. 

Let us designate the mass function at scales v 3> 1, i.e, 
in the classical regime, as F = (87rG /3)pi? 3 , whereas 
for v Ri a^ 1 (in the quantum geometry regime) we use 
F e ff given by Eq. (3.20). Then, the mass loss, AF/F, is 



provided by the following expression: 

AF F eS = p 

F F Pc 



Using Eq. (3.19) in Eq. (4.14), we can estimate the time 



F 2 



variation of the mass function as 



(4.10) 



F 



cfr 



-8ttG q HR a p 1-3 



8nG R 3 



Pc 

3 J \ p c 



8nG 



P 1 



(4.15) 

The classical Raychaudhuri equation, 2 = — 47r ^' (p + 
3p), in effective theory can be modified from Hamilton's 
equation as follows: 



4ttG q 



-2p 1 



+ 6p 1 



>1_ 

Pa 



(4.16) 

So that, by rewriting the modified Raychaudhuri equa- 

■ ttG 
3 



tion as § = — 4 , (peff+3p c ff ) ; we can define the effective 



pressure as 



This equation can be interpreted as follows: If p <C per, 
i.e. at the classical limit, then AF « 0. However, when 
the energy density p — )■ p cr at the critical value (bouncing 
point) in the semiclassical limit, we have AF/F — » 1; 
this means that when the collapsing cloud approaches 
the critical energy density, an outward flux of energy will 
start near the bounce. 

It is worthwhile to mention that, within an inverse 
triad correction of collapsing system (with a scalar field 
[TU] . or a tachyon field [13], as matter source), the 
(quantum) modified energy density decreases as collapse 
evolves. Whence, as the collapsing cloud approaches the 
center (with a vanishing scale factor, where the classical 
singularity is located) the energy density reaches its min- 
imum value, whereas the mass loss tends to one. In the 
model herein, as the gravitational collapse proceeds in 
the quantum regime, the energy density increases and 
reaches its maximum value p CT at the bounce (with a fi- 
nite non-zero volume). In this process, however, the mass 
loss tends to unit earlier, and before that the collapse 
reaches to the center of star. 

Let us assume that the energy density flux is meas- 
ured locally by an observer with a four-velocity vector 
Then, the energy flux and the radiation energy 
density are measured in this local frame and given by 
cr = Tfj, v ^^ v . Furthermore, the total luminosity for a ra- 
dially moving observer with radial velocity 1? = £ rv = 
at the radius r v , is given by L(v) — Anr%(J [35], so that: as expected from Eq. ( |2.3| . Eq. ( |4.18[ ) shows that in last 



p eff := p 1 - 3 



Pc 



(4.17) 



Inserting Eq. (4.17) in Eq. (4.15), the time derivative of 
the mass function can also be expressed as 



F eff = 8irG R 



/8ttG 
l~3" 



Peg PcS- 



The classical limit of Eq. (4.18) is 



F = 8nG QP R 3 f 8 ^ -V - — 



(4.18) 



p = -8ttG pR z R, (4.19) 



L(v) = - 



dM(v) 



(7 + 19) 2 dv 



(4.11) 



where 7 = (l+i? 2 -2M(v)/rv) _1 . As long as dM/dv < 0, 
the total luminosity of the energy flux is positive; this 
indicates that there exist an energy flux radiated away 
from interior spacetime and reaching the distant ob- 
server. Moreover, from = — 1 and £ e = £ v = 0, 
we have 



stage of the gravitational collapse, the mass function de- 
creases when the effective pressure becomes negative. In 
contrast, the classical system shows a continuous increase 
of the mass function, with an increasing positive pressure, 
pointing to a black hole end stat e. In figure [4] we have 
a graphical representation of Eq. (4.18 ) during evolution 



dv 



= e 



1 

7 + i9 



(4.12) 



Now, using Eqs. ( 4.11[ ) and ( 4.12[ ), and substituting 
M(v) with M = F(t,n)/2 at the boundary E we ob- 
tain, 



L(v) = - 




(4.13) 



For an observer at rest (1) = 0) and infinitely distant 
(r v — > 00), the total luminosity of the energy flux , L^, 
can be obtained by taking the limit of Eq. (4.13): 



£00 (*) = - 




(4.14) 



of the collapse. Therein, we have that the effective mass 
function's time derivative, being initially positive (like its 
classical counterpart), starts to decrease until it becomes 
negative and then vanishes at the bounce. At this stage, 
as Lqo ~ — F c s, the luminosity positiveness is related to 
the mass loss occuring near the bounce. The existence 
of a negative F e g, with some substantial variation in the 
vicinity of the bounce, points to the existence of an en- 
ergy flux radiated away from the interior spacetime and 
reaching the distant observer [25] . 

Fig ure [S"| presents the behaviour of the effective pres- 
sure (4.17) conveniently scaled with the critical density 
p cr . Therein we should stress that the effective pres- 
sure becomes super negative near the bounce and takes 
the value p c ff(Pcr) = — 2p cr at the bounce. This verifies 
the general understanding derived from homogeneous iso- 
tropic models indicating that, the singularity resolution 
is associated with the violation of energy conditions. This 
suggests that quantum gravity provides a repulsive force 
at very short distances [25] . 
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4 
a 



Figure 4: Behaviour of the time derivative of the mass func- 
tion, F, in the effective quantum dynamics regime (from Eq. 
(4.18l) for the value of parameters Go = ciigw = 1, and 



7T0 = 10 000. The full line is for effective dynamics, whereas 
the dashed line is for the classical counterpart. 



2^ 



1 - 



Peff 
Per 



-1 



4 
a 



Figure 5: Behaviour of the effective pressure (from Eq. ( 4.17 1) 
for the value of parameters Go = cu g ht = 1, and — 10 000. 
At the bounce we have jJ e ff(p C r) = —2p CI . The full line is for 
effective dynamics, whereas the dashed line is for the classical 
counterpart. 



B. Non-singular black hole formation 

If the initial mass of the collapsing star is bigger than 
the threshold mass m*, given by Eq. (4.9), then a black 



hole would form at the collapse final stage. In this sec- 
tion we will analyse a possible prediction for the exterior 
geometry of the collapsing system in this case. 

The total mass measured by an asymptotic observer 
is given by m oxt = tum + m$, where ttlm is the total 
mass in the generalized Vaidya region, and — j pdV 
is the interior mass related the scalar field <f>. Since the 
matter related tum is not specified in the exterior Vaidya 
geometry in our model, we just focus on a qualitative 
analysis of behaviour of the horizon close to the matter 
shells. 

From the matching conditions \2.11 )-(2.14), we can get 
the information regarding the behaviour of trapping ho- 
rizon in the exterior region. Indeed, when the relation 



2M(v,r v )G = r v is satisfied at the boundary, trapped 
surfaces will form in the exterior region close to the mat- 
ter shells. In classical geometry, the boundary function, 
J- = (1 — 2M(v)Go/r v ), becomes negative in the trapped 
region and vanishes on the apparent horizon. Therefore, 
the equation for event horizon is given at the boundary 
of the collapsing body by — 0. Nevertheless, in semi- 
classical regime, the boundary function is expected to be 
modified by employing the matching conditions due to 
the fact that the interior spacetime was modified by the 
quantum gravity effects. 



Using the conditions (2.11) and (2.13), the classical 



function J- can be written as 



71 



1 



F 
R ' 



(4.20) 



where we have used 2M(v,r v )Go/r v — F(t)/R(t) at the 
boundary surface £ with r — rv In semiclassical regime, 
the mass function is modified as Eq. ( 3.19[ ) , whereby the 



boundary function (4.20) is modified as 



F 
R 



1 - 



F 2 



(4.21) 



at the boundary shell r&. Eq. (4.21) shows that the 



quantum gravity effects leads to a modification of the 
boundary function by a cubic term F . 

By substituting the classical mass f unctio n with F = 
(8irG /3)pR 3 , we can rewrite the Eq. (4.21) as 



^(*) = i-^(i-Jr 



(4.22) 



where A = (47rG /3)7r^r^, and B = 7r^rf /(2p cr ) are con- 
stants. Eq. ( 4.22[ ) represents a non-singular, exotic black 
hole geometry. Notice that, this function for large values 
of R tends to the classical limit: 



F{R) = 1 



.4 
R 1 



(4.23) 



which represents a classical singular black hole geometry 
|12j . In this case, the apparent horizon could form 
when J- in Eq. (4.23) vanishes; equation of this hori- 



zon which intersects with the matching surface is given 
by a(t) — A 1 / 4 /r b . Furthermore, as collapse tends to the 
singularity as R — > 0, the boundary function diverges 
as T — » — oo; the classical singularity is covered by the 
black hole horizon. In addition, as we expected, in the 
presence of a nonzero matter pressure (of the massless 
scalar field) at the boundary, the (homogeneous) interior 
spacetime can not be matched with an empty (inhomo- 
geneous) Schwarzschild exterior pTj . 



V. CONCLUSIONS AND DISCUSSION 

We considered a spherically symmetric and homogen- 
eous spacetime for a collapsing system, with a massless 
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scalar field as matter content. The homogeneous interior 
is matched with an exterior Vaidya geometry. 

We subsequently studied the interior spacetime within 
the effective theory of loop quantum gravity provided 
by a holonomy correction to the dynamics of the system 
|18j . It was shown that loop quantum effects removes 
the classical singularity that arises at the end state of 
the collapse, and replaces it by a quantum bounce. Fur- 
thermore, we investigated the evolution of the trapped 
surfaces emerging from the semiclassical interior space- 
time. The physical modifications related to the semi- 
classical regime provided three cases for the trapped sur- 
faces formation, depending on the initial conditions of 
the collapsing star. In particular, our solutions showed 
that, if the initial mass of the collapsing cloud is less 
than a threshold mass, no horizon forms during the col- 
lapse, whereas for the mass equal and larger than the 
threshold mass, one and two horizons form, respectively. 
It is worthy to mention that, this scenario is qualitatively 
similar to the model previously predicted from an inverse 
triad correction [TT] . 

The interior semiclassical effects is carried out to the 
exterior geometry by the matching conditions on the 
boundary of two regions. Hence, an effective geometry 
emerged for the exterior metric which describes the phys- 
ical consequences for the late-time evolution of the col- 
lapse. In the case in which no horizon forms, we have 
showed that, as the collapse evolves, the energy dens- 
ity and mass function increase towards the maxima p cr 
and F CI , respectively. However, this energy growth is 
accompanied by the effective mass loss, induced by the 
quantum geometry regime. Therefore, in the late-time 
stages of the collapse, it gives rise to an outward flux of 
energy from the interior quantum spacetime to the dis- 
tant observer. Similar results was discussed for a scalar 
field collapse, where an inverse triad modification was 
employed |10j . However, the mass loss obtained therein, 
was characterized by a reduction of the energy density 
and mass function towards the center of the star. In the 
other case in which one or two horizon form, the exterior 
geometry obtained predicts an exotic black hole forma- 
tion that is different than the Schwarzschild one. 

It has been proposed that gamma-ray bursts may oc- 
cur at the final state of the collapse [27 . In our model 
herein, although the effective theory is not accurate in a 
deep Planck regime, there might be some stages of the 
gravitational collapse where the effective quantum geo- 
metry becomes relevant. Therefore, from such a scenario 
it could be observed some of those outward flux of energy 
from gamma-ray bursters. 

Modification to the effective mass function through 
the quantum geometry effect at the boundary of two re- 



gions, results in an effective exterior Vaidya geometry. 
More precisely, from Eq. (3.201 we have that F e g — 
F(G eS /G ) . Therefore, using Eq. ((2H3]), we get F cS = 
2MG cS , so that 



M eff G - MG eS . 



(5.1) 

Consequently, the Vaidya metric (2.10) takes the im- 
proved form: 



g+dx^dx" = - (1 - 2M(v)G cS /r v ) dv 2 



2dvdr v 



:dn 2 . 



(5.2) 



The improved Vaidya solution (5.2 1 shows that the 
effective quantum (interior) spacetime is manifest 
through a modification to the Newton constant Go as 
Go (1 — pi Per)- Similar scenario was presented in recent 
works [28-30]. Therein, it was shown that, in situations 
where the interior spacetime undergoes a transition from 
positive to negative pressures, a quantum improvement 
to the Vaidya outgoing solution should be implemented. 
Nevertheless, this improved outgoing Vaidya solution, 
results in a modification of the Newton constant as G (r v ) 
being a scale dependent function [2S]. The redefinition 
of the Newton constant in this context was introduced to 
model quantum effects and avoid the problem of unbound 
back scattered radiation [2"8] . 

The qualitative picture that emerges from our toy 
model was influenced by a homogeneous interior space- 
time, however, in a realistic collapsing system, one has 
to employ a more general (inhomogeneous) setting. Non- 
etheless, in the full dynamics of states in LQG, it seems 
that the 'energy' in the inhomogeneities is sufficiently 
small even in the Planck era, where the inhomogeneous 
modes will be affected by the quantum geometry effects 
of the dominant homogeneous modes |31j . In this case, 
the inhomogeneous modes would be too weak for their 
own quantum geometry effects to be important. There- 
fore, the effective theory, as an approximation of the im- 
proved dynamics LQC, can be a good phenomenological 
approach to study the gravitational collapse. 
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